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Abstract 

The evolution of a dispersal trait is a classical question in evolutionary ecology, which has been 
widely studied with several mathematical models. The main question is to define the fittest dispersal 
rate for a population in a bounded domain, and, more recently, for traveling waves in the full space. 


In the present study, we reformulate the problem in the context of adaptive evolution. We consider 
a population structured by space and a genetic trait acting directly on the dispersal (diffusion) rate 
under the effect of rare mutations on the genetic trait. We show that, as in simpler models, in the 
limit of vanishing mutations, the population concentrates on a single trait associated to the lowest 
dispersal rate. We also explain how to compute the evolution speed towards this evolutionary stable 
distribution. 

The mathematical interest stems from the asymptotic analysis which requires a completely different 
treatment for each variable. For the space variable, the ellipticity leads to the use the maximum 
principle and Sobolev-type regularity results. For the trait variable, the concentration to a Dirac 
mass requires a different treatment. This is based on the WKB method and viscosity solutions 
leading to an effective Hamiltonian (effective fitness of the population) and a constrained Hamilton- 
Jacobi equation. 


Key words: Dispersal evolution; Nonlocal pde; Constrained Hamilton-Jacobi equation; Effective 
fitness; 
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1 Evolution of dispersion 

Evolutionary dynamics of a structured population. There are several well established math¬ 
ematical formalisms to describe evolution. Game theory is widely used since the seminal paper [29]: 
see also [23] . Dynamical systems are also employed to describe the possible invasion of a population 
be a mutant, and to characterize Evolutionary Stable Strategies and other mathematical concepts; 
see, for instance, mm- Stochastic individual based models are often used to describe the evolution 
of individuals undergoing birth, death and mutations. Relations with other approaches can also been 
made at the large population limit, see El 113 and the references therein. 
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The formalism we use here is still different and relies on a population structured by a phenotypical trait 
and competing for a limited resource. This approach was initiated and has been widely studied in [El- 
Several other versions use the formalism of a population structured by a trait, undergoing mutations 
and competition [Ml 021 [251 [28]. All these papers, however, consider only proliferative advantage. The 
population with the highest birth rate or best competition ability survives while other traits undergo 
extinction. Mathematically this is represented by a limiting process where the population number 
density, denoted here by to, takes the form of a weighted Dirac concentrated at the fittest trait 9. 

The extension to models, where the phenotypical trait is combined with another structuring variable, 
usually space, is more recent and leads to considerable mathematical difficulties; see m m EziGni na 
[10] . This is mainly due to the fact that for the trait variable 9 the solutions concentrate as described 
above while in the space variable solutions remain bounded. 

Another motivation for considering the model in this note is to study the selection of the fittest 
individuals without a proliferative advantage. In this context, the reproduction rate might be com¬ 
pensated by another advantage. This gives rise to the question of defining an “effective fitness.” The 
gradient of the effective fitness determines the direction of trait evolution and its maximum defines 
the evolutionary stable strategy. 

A particularly interesting example, both mathematically and biologically, in this directions is the 
selection of a dispersal rate, which we describe next in the context of a continuous dispersal trait. 
When only two species are represented by their number densities toi(x), TO2(x) and are competing for 
the same resource (carrying capacity) I\, the question is to know which of the two dispersal rates D \, 
Z ?2 is prefered in the competiton. Then, the model is 

— DiA x m = m[K(x) — toi + TO 2 ), ie!l, i = 1,2. 

Is it better to ‘move’ faster or slower? In other words, is it favorable to have D\ > D2 or the contrary? 


The model. Here we assume that all dispersal rates are possible, and include mutations. We ask 
the question of the rare mutation limit of the steady state version for a still simple model of evolution 
of dispersal in a population. The main modeling assumptions are: (i) all individuals wear a phenotype 
characterized by a parameter 9, which induces a dispersal rate D{9 ), (ii) a Fisher-type Lotka-Volterra 
growth/death rate with a space dependent carrying capacity K and limitation by the total population 
whatever the trait is, and (iii) rare mutations acting on the genetic variable and modeled by a diffusion 
with covariance \f2e\ we refer to }12| for a derivation of this type of equations from individual based 
stochastic models. 

More precisely, we study the asymptotic behavior, as e —> 0, of the density n £ = n e (x, 9) 

(0,oo), with Slci^a smooth domain, of the nonlocal and nonlinear problem 


-D(9)A x n e - e 2 A e n £ = n e (K{x) - p £ {x)) in Qxl, 

< f 1 (!) 

p £ {x) = / n e (x, 9)d9, 

Jo 


with Neumann boundary condition on <9fl and 1— periodicity in 9 , that is, if u is the external normal 
vector to D, 


d 

— — n £ = 0 on dfl X M and n £ is 1 — periodic in 9. 
ov 


( 2 ) 


We have chosen periodic boundary conditions in 9 to simplify some technical aspects concerning a 
priori estimates. 
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As far as the carrying capacity K and the dispersion rate D are concerned, we assume 

K G C' 1 (0), there exists K m > 0 such that K(x) > K m > 0, K is not constant. (3) 

and 

f DgC 1 (M; (0, oo)) is 1-periodic and 

I there exists a unique local minimizer 9 m G [0,1) such that D > D(9 m ) := D m > 0; 

we note Q is used to assert that the effective Hamiltonian also has a minimum at 9 m . 

Finally, for technical reasons, we also assume that 

H is convex. (5) 

It follows from, for example, mm, that, given © and Q, the problem Q and ([2]) admits a strictly 
positive solution n £ : H x M —> (0, oo). 


Formal derivation of the mathematical result. We proceed now formally to explain what 
happens in the limit e —> 0 and, hence, motivate the statement of the results. As it is often the case 
with problems where it is expected to see concentration in the limit, we make the exponential change 
of variables 

n £ = exp (u £ / e), (6) 

which leads to 


m 




D(d) 




K(x) — p e {x) in fixl, 


(7) 


with 


Ur is 1—periodic in 9 and — u £ = 0 on dQ X M. (8) 

ov 

It is clear from 0 that, if the u £ s have, as e —X 0, a limit u , it must be independent of x. and, hence, 
it is natural to expect the expansion 


u £ (x, 9) = u{9) + e lnA7(x, 9) + 0(e). (9) 

Assuming that, as e —> 0, the p £ s converge to some p , a formal computation suggests that A f is the 
positive eigenfunction of 


f-D(9)A x N = N (K~p) +A r H (9, p(-)) in H, 

\ = 0 on dVL, 

with eigenvalue H[9,p(-) ), and that u solves the constrained Hamilton-Jacobi equation 

-\V e u\ 2 = -H(9, p(-)) in M, 


< maxflei u{9) = 0, (11) 

u is 1— periodic. 

The constraint on the maxu becomes evident from the facts that, as it turns out, the p £ s are bounded 
uniformly in £ and the equality 


p e (x) = f n £ (x, 9)d9 = f exp (u £ (x,9)/e)d9, (12) 

Jo Jo 

which also suggests that, as e —> 0, the n e ’s behave like a Dirac mass with weight p. 
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The mathematical result. To state the result we recall that 5 denotes the Dirac mass at 0 and 
we introduce the nonlinear Fisher-type stationary problem 

f -D m A x N m = N m (K - N m ) on H 

l lb N ™ = 0 on <9D, 

which, in view of ([3]) and 0. admits a positive solution N m > 0 (see, for example, IE!-) 

We have: 

Theorem 1.1 Assume ([3]), Q, and ©■ Then, as e ^ 0 and in the sense of distributions, 

n e —¥ N m (x)d(9 — 6 m ) and p e ^N m . 

Moreover, as e ^ 0 and uniformly in x and 6, eln n £ —> u, where u is the unique 1—periodic solution 
to with p — N m , such that m9.x u — — 0. 


Biological interpretation. The conclusions of Theorem 0 can be thought as a justification of 
the fact that the population selects the “slowest” individuals in accordance with several previous 
observations on the evolution of dispersal. In this respect, the eigenvalue —H(Q,p{-)\ defines the 
fitness of individuals depending upon their trait. This fact can be stated using the canonical equation 
(31), which is formally derived in Section[6j In the words of adaptive dynamics, our result characterizes 
the unique Evolutionary Stable Distribution (or Strategy), |16l 25] , 


That mutants with lower dispersal rates can invade a population, that is the characterization of 6 m by 
property © , is known from the first mathematical studies [221 H8j. However, these papers use time 
scale separation, heuristically assuming a mutant appears, with ‘small mutation’, and compete with 
the resident population. Our approach here is more intrinsic since we consider structured populations 
competing for resources and undergoing mutations. Surprisingly when set in the full space where the 
problem is characterized in terms of traveling waves, the opposite effect is observed, that is mutants 
with higher dispersal rates are selected giving rise to accelerating waves, 0 El EH HUS]. For two 
competing populations, the combined effect of dispersal and a drift is studied in m- The analysis of 
dispersal evolution also gave rise to the notion of ideal free distribution mm- 


The question of dispersal evolution is a classical and important topic in evolutionary biology. The 
reader can consult |3B] for a survey of the many related issues, to m for the case with patches and 
demographic stochasticity, to [24] for the case of trajectories with jumps (nonlocal operators) and to 
[H] for other biological references about accelerating fronts. A formalism using Fokker-Planck equation 
is used in [35]. Also, let us mention that a remarkable qualitative aspect fo spatial sorting is that in 
the full space, the largest dispersion rate is selected mum- Finally, another mathematical approach 
to the concentration effect, stated in Theorem 0 can be found in [26]. 


Organization of the paper. In Section [2] we prove some uniform in e estimates for the p e ’s that 
are then used in Section [3] where we derive the effective Hamiltonian, that is the eigenvalue problem 
of (10). In Section [4] we introduce the constrained Hamilton-Jacobi equation to conclude the proof 
of Theorem |1.1| and Theorem |4.1| In Section [5] we prove the two technical lemmata that were used 
in Section [4] Finally, In Section [6j we provide some perspectives about the problem, namely a more 
precise asymptotic expansion and the parabolic case, as well as numerical examples for the evolution 
driven by the parabolic equation. 
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2 Estimates on p £ 


We state and prove here some, uniform in e, estimates for the p e ' s, which are fundamental for the 
analysis in the rest of the paper; here |fl| is the measure of 0 and C(A) denotes a constant that 
depends on A. 


Lemma 2.1 Assume ([3]) and Q. There exist positive independent of e constants C\ = C\{K,D,VL), 
C 2 = C 2 (K, D, Cl) and C 3 = C 3 (|f 2 |, min/\) such that 

(i) 0 < p £ < Ci, 

< (ii) sup ee(0il) \\p e \\w 2 ’P(n) < C 2 for all p £ [l,oo), and (14) 

,(«) f Q pe{x)dx > C 3 , 

and, along subsequences £ —> 0 , the p £ ’s converge uniformly in O to p £ C(Q). 


Proof. We first observe that p £ trivially satisfies the Neumann condition 

d 


du 


p £ = 0 on dQ. 


After dividing 0 by D(8), integrating in 8 and using the periodicity in 9, we find 


1 


1 


-A x p e -£ n £ (x,9)Ag—d8 = / -n £ d9(K-p £ ), 


D 


D 


and, hence, for some constant C , which only depends on D and K , we have 


A x p £ + 


1 


\D\ 


Pi < c Pe . 


(15) 


Then ([l4|)( i) follows from the strong maximum principle, while the W 2,p -estimates are a consequence 
of the classical elliptic regularity theory. 

The lower bound (|14[)(iii) comes from integrating ([!]) in x and 9. Indeed, in view of the assumed 
periodicity, we find 


max K 


p £ (x)dx > I p s {x)K{x)dx = f p £ {xfdx > tttt ( / p e (x)dx) . 
! Jn Jn |“| \Jn ) 


The last claim is an immediate consequence of the a priori estimates and the usual Sobolev imbedding 
theorems. 


□ 


3 The effective Hamiltonian 

For 8 £ [0,1] and p £ L°° (0; [0, 00 )) we consider the eigenfunction AT = M(-', 9, p(-)) £ iT(O) and the 
eigenvalue H = H(6,p(-)\ of 

f-D{9)A x Af = AT (K - p)+M H in 0, 

|j^AA = Oon<9f2 and f Q N{x] 9, p(-)) 2 dx = 1 . 
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Note that in view of ([3]) and the regularity of ff, the existence of the pair (jV(-; 9, p(-)), H(Q, p(-))) 
follows from, for example, mm- 

The next lemma provides some important estimates and information about 77(0,/?(•)). In the state¬ 
ment we use the notation Km ■= max^ K. 

Lemma 3.1 Assume Q and Q. Then 

(*) - K m < H(6,p(-)) and (u) H{9, p(-)) < t^j p(x)dx, (17) 

the maps 6 —> H{9, p(-)) and 6 —> D{9) have the same monotonicity properties for all p , (18) 

and, in particular, 9 m is the unique minimum of H(9, p) in [0,1] for all p. 

Proof. Multiplying the equation by Af and integrating in x gives 

0 <D{9) [ |VA7| 2 = [ N 2 [K-p + H]dx< [ N 2 dx[K M + H], 

Jq Jq Jq 

and, thus, pT|)(i) holds. 

The upper bound © (ii) follows from the positivity of K , since, after dividing the equation by A f 
and integrating by parts, we find 


f -— 1 ro - dx = [ [K — p + H]dx< 0. 

In A7 2 Jn 1 


For (18), we differentiate in 9 the equation in (16) to find 


- D'(9)A x Af - D(9)A x Af e = N e (K(x) - p{x) + H) + N H g , 


(19) 


where Afg and Hg denote derivatives with respect to 9, we multiply by M and integrate by parts using 
the boundary condition to get 


D\6) [ |V X A r\ 2 + D(9) [y x M.V x Me= (NN e [K-p\ + H e {e-p{-)). 
Jn J J 


( 20 ) 


Next we use the fact that the L 2 — normalization of M yields j^AfgMdx = 0 and after multiplying the 
equation in (16) by Afg, integrating by parts and subtracting the result from (20) we find 


D'{9) [ \V x Af\ 2 = H e {9-p{-)). 
Jn 


Since we have assumed in ([3]) that K is not constant, V X A7 does not vanish and the result follows. 

To conclude, we observe that Q and (18) yield that H as the same monotonicity, in 9, as D and, 
thus, has a unique local minimum at 9 m . 

□ 
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4 The constrained Hamilton-Jacobi equation 

We prove here the generalized Gaussian-type convergence asserted in Theorem [13 derive the con¬ 
strained Hamilton-Jacobi equation © and state some more properties. For the benefit of the reader 
we restate these assertions as a separate theorem below. 

Theorem 4.1 The family u £ is uniformly in £ Lipschitz continuous and converges, uniformly in x and 
9, to u, which is independent of x and satisfies, in the viscosity sense, the constrained Hamilton-Jacobi 
equation ©• Moreover, 

min H(9,p(-)) = 0 = H(9 m ,p(-)) = H 0 (6 m ,p( •)). (21) 

6 


Since —H represents the fitness, it turns out that © characterizes 6(6 — 9 m ) as the Evolutionary 
Stable Distribution (or 9 m as the Evolutionary Stable Strategy). See 


It follows from (21) that both H and its derivative vanish at 9 m . As a result the viscosity solution u 
of © also vanishes at 9 m as u(6) = O((0 — 9 m ) 3 ^ 2 ), and this makes the connection with the result 
of 


The proof is a consequence of the next two lemmata which are proved later in the paper. 
Lemma 4.2 (Bounds on u £ ) There exists an independent of e constant C such that 


ma xu e (x,9)dx < Ce. 


( 22 ) 


Lemma 4.3 (Lipschitz estimates) There exist an independent of e constant C such that 

-rrl V x u £ \ 2 + I Vgu £ \ 2 < C and max u £ (x,8)<Ce. 

E 1 a;ef2,0<6Kl 

Moreover, the u £ ’s converge, along a subsequence £ —> 0 and uniformly in x and 9, to a Lipschitz and 
1-periodic function u : M —> M such that maxo<6><i u(9) = 0. 


We continue with the proofs of Theorem 4.1 and Theorem 0 

Proof of Theorem |4.1| . The fact that any limit u of the u £ s satisfies © is a standard application 
of the so-called perturbed test function method (see j[[9j) and we do not repeat the argument. 

It follows from (11) that H ( 9 , p(-)) > 0, while, at any maximum point 9 of u, we must have H(9, p (•)) < 
0, and, hence, H(9, p(-) ) = 0, and, in view of Lemma 


3.1 


min H ( 9 , p(-)) = 0 = H(8 m , p(-)). 

u 


As a result the only possible maximum point of any solution of © must be at at 9 m , which implies 
that the equation has a unique solution. 


Also the knowledge of 9 m determines uniquely the limit p(x) 
the full family ( p £ ,u £ ) converges. 


N m (x ), from equation (13), and, thus, 


□ 
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Proof of Theorem |1.1| . The statement in terms of n £ is an immediate consequence of Theorem |4.1 
Because u e achieves a unique maximum at 9 m , from the Laplace formula for n e written as ©> , we 
conclude that the n E (x , 0)’s converge weakly in the sense of measures to p(x)5(9 — 9 m ), with p{x) the 
limit of p £ (see Section [2]) . 

Next, integrating equation Q in 9 we conclude that 


-A 


D(9)n e d9 = p £ (x)(K(x) - p £ (x)). 


Passing to the limit e —^ 0, and taking into account that n is a Dirac mass at 9 m , we find that p = N„ 


because they both satisfy the equation (13). 


□ 


5 The proofs 

We begin with the proof of Lemma |4.2| 

Proof of Lemma 4.2, Integrating (J7J), we find 


-oil \'\7 x u £ \ 2 dxd9 + / / \\7 g u £ \ 2 dxd9 < C. 


(23) 


The claim for the maximum bound follows from the L°°— estimate on p £ . 

For each iGll, let M e (x) := ma xqu £ (x,9) and choose 9 e {x) such that M e (x) = u e (x,9 e (x)). Then it 
follows from (|14[)(i) that 


M, 


(x)/e f 1 e (u e { x ,e)-u e (x,0 e ))/s d6 = Pe ^ < C 

Jo 


Inserting in (24) the estimate 


r e e( x ) /-i 

u e (x,9) - u e (x,9 e (x)) = / d e u £ (x,9')d9'> -1- \d g u £ (x, 9')\ 2 d9', 

Je Jo 


we get 


and, hence, 


--[1+ / \d g u £ (x, 9')\ 2 d9'j 


e M e (x)/e & £ j o 


<C, 


Using (23), we obtain 


- -[1 + f 1 \d e u £ (x,9')\ 2 d9'\ < In C. 

£ £ Jo 


M e {x)dx < Ce , 


(24) 


and this concludes the proof. 


Next we discuss the proof of Lemma 4.3 


□ 


Proof of Lemma 4.3|. We first assume the Lipschitz bound and prove the rest of the claims. 














Let M £ := ma x xG qM £ (x), with M e (x) as in the proof of the previous lemma. It is immediate from 
p2] ) and the fact that \V x u £ \ 2 < eC that, for some C > 0, 

M e < eC. 


Next we show that lim inf^o M e > 0, which, in turn, yields that lim^o-Me = 0. 
Let ye 6 be such that M e := M e (y e ) = u £ (y e , 0 e (y £ )) and write 


p e ( X )=e M '/ e f 1 e (u e (x,e)-u e (y e A))/s de 

Jo 


Combining the lower bound on p £ in Lemma 2.1 and the (Lipschitz) estimate |V x ri £ | 2 < eC, we get 

C 3 < [ p £ (x)) <\n\e M ^ e e c [\( Ueiye ’ e) - u ^ ysA) ) /£ d9 <e M ^ E e c \n\, 

in Jo 

and, thus, M e > —Ce. 

Now we turn to the proof of the Lipschitz bounds, which is an appropriate modification of the classical 
Bernstein estimates to take into account the different scales. We note and prove Lemma 4.3 Note 
that the convexity assumption on 17 is used solely in this proof. 

We begin by writing the equations satisfied by \D x u e \ 2 and \dou £ \ 2 which we obtain by differentiating 
0 in x and 9 and multiplying by D x u e and d g u £ . We have: 


— ^p-A x \D xUe \ 2 — eA g \D x u £ \ 2 + 2 : ^|Z7 2 3 ,i( £ | 2 + 2£|L» 2 £) rt £ | 2 

-2*JJp-D x u £ .D x \D x u £ \ 2 - 2d g u £ .d g \D x u £ \ 2 = 2D X (K - p £ ).Du £ , 


(25) 


and 


12 "iel 2 


A x \d e u £ \ 2 - eA g \d g u £ \ 2 + 2JJJp-\D 2 9 u £ \ 2 + 2e\Dj d u, 

-2JJJp-D x u £ .D x \d e u £ \ 2 - 2d g u £ .d g \d g u £ \ 2 = 2 I JJp-A x u £ d g u £ + A^J\D x u £ \ 2 d g u £ . 

Let q = Dx £ 2 E - + \d g u £ \ 2 and compute 

- JJJp-A x q - eA e q + 2^-\Dl x u e \ 2 + 1+D £ [e) \D 2 x6 u £ \ 2 + 2e\Dj e u £ \ 2 

-2 JJJp-D x u £ .D x q - 2d g u £ .d g q = 4 D X {K - p £ ).Du £ + 2 ^M A x u £ d g u £ + A^JP-\D x u £ \ 2 d g u £ . 


(26) 


Assume that (x,9) is a maximum point of q. Because of the convexity assumption, we have J^q < 0 
on the boundary (see W) and thus x £ ^ <912 X [0,1]. 

Therefore, at this point (x,9), we have 

^\D 2 xx u £ \ 2 + e\D 2 ee u £ \ 2 < XD x (K - p e ).Du e + ^-A x u £ d e u £ + 2^ r \D x u £ \ 2 d g u £ 

— ^D x (K — p £ ).Du £ + 5^(A x u £ ) 2 + | {d g u £ ) 2 + \D x u £ \~d g u £ 
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and we choose 5 small enough so that we can absorb the term 5^(A x u £ ) 2 in the left hand side. 
Since there is a constant C±(d, D) such that 


— A x u £ + eA e u £ 
£ 


< Ci^t\D xx U £ \ 2 + Cie\DggU £ \ 2 , 


we conclude (using the equation) that, for some C*,{d,D,5) > 0, 

12 i2 


+ + K _ pt 


<c 2 


1 2 

~Pe)‘D u £ T e{deu £ ) -1- dQii £ 


It follows that there exists some positive constant C such that 

q 2 < C [l + \\Du e \ +E 2 (d e u £ ) 2 + 2 \D x u £ \ 2 \d e u £ \] 

< c [l + y/q + e 2 q + eq 3/2 ] . 

From this we conclude that q is bounded and the Lipschitz continuity statement is proved. 


□ 


6 Conclusion and perspectives 

Conclusion We have studied a steady state model describing the Evolutionary Stable Distribution 
for a simple model of dispersal evolution. The novelty is that mutations acting on the continuous 
dispersal trait (the diffusion rate) are modeled thanks to a Laplacian operator, and this replaces the 
standard ‘invasion of a favorable mutant’ in the usual time scale separation approach. When the mu¬ 
tation rate is small, we have shown that the minimal dispersal is achieved, in accordance with previous 
analyses. This is mathematically stated as the limit to a Dirac mass which selects the minimum of 
the diffusion coefficient in the equation. Technical difficulties rely on a priori estimates in order to 
make the approach rigorous and establish the constrained Hamilton-Jacobi equation which defines the 
potential in the Gaussian-like concentration. 

We indicate two possible extensions of our results. The first concerns the way to make more precise 
the convergence result of Theorem 0 The second is about the time evolution problem. 


A more precise convergence result. The question we address here is whether it is possible to 
make more precise the convergence n £ —> N m (x)5(8 — 9 m ) in the weak sense of measures stated in 
Theorem 11.11 


The gradient bound in Lemma 4.3 implies that || u e (x,9) — u £ (0, 0)||oo < C\\'Vu £ (x,0\\ oo < Ce, and, 
therefore along subsequences, the ’ s converge in L°°— w*. To prove more about the 

corrector, it necessary to have an estimate for u £ (O,0) — u(9), which is not, however, available from 
what we have here. 


Another approach is to introduce, instead of u(9), the eigenvalue problem 


-£ 2 A e W £ = W e [-H(d, p £ ) + A e ], W £ is 1-periodic, W £ > 0. 
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The change of unknown w £ = e In W £ gives the equation 

-eA g w £ - \D g w £ \ 2 = —H(0, p £ ) + A e . 

Standard gradient estimates yield that \Dqw £ \ is bounded independently of e, and, from the equation, 
we see that the eAqw £ s are also bounded independently of e. It follows that, as e —» 0, A e —> 0 and, 
after an appropriate normalization, w £ —> u uniformly. 

As a result we can factorize the solution of 0 as 

n e (x,9) = v £ (x,9)e w * (e) / £ , 


and we claim that, for some other factor p, 


v £ (x, 9) —> pJ\f(x, 9), L°°-w-k, 


(27) 


a statement which is more precise than that of Theorem 0 
To see this, we write the equation for v £ as 


-A x v £ - e 2 Aqv £ - 2 eDgWe ■ D e v £ = v £ [K - p £ + H(9, p e ) - X £ j. 


To obtain bounds on v £ we notice that p £ (x) = v £ W £ dO is bounded (from above and below), and, 
since W £ —>• 5(9 — 9 m ), we conclude that, for some 0 £ and (L near 9 m and independent of e constants 
C and c, v £ (x,0 £ ) < C and v £ (x,9 e ) > c. It then follows from standard arguments that the u £ ’s are 
bounded in L°°. 


The limiting equation is the eigenfunction problem (16) and positive solutions are all proportional to 
A f(x,9), which gives the statement. 


The parabolic problem. Our approach is mainly motivated by the dynamics of the evolution of 
dispersal. The steady states are the Evolutionary Stable Distribution and are obtained as the long 
time distribution of competing populations |25j. This leads to the study of the time dependent problem 


en £: t(x , 9 , t ) — D{9)A x n e — e 2 Ag = n £ (K(x) — p e (x, t )) in 0 x M x (0, oo), 
p £ (x,t) = f* n £ (x, 9, t)d9, 


(28) 


with the Neuman boundary conditions on dfl and 1—periodicity in 9. 


For this problem there are two limits of interest, namely e —> 0 and t —> oo. So far, we have studied 
the limit t —> oo, reaching a steady state ([!]) of (28), and then considered the limit e —> 0. 


Reversing the order, we need to study first what happens as e —> 0. In this case, we expect that, in 
the weak sense of measures, 


n e (x,9,t) —>• N(x,t)5(9 — 9(t)) and p £ —> N(x, t), 
where at least formally, the weight N(x, t ) is defined, for each t, by the stationary Fisher/KPP equation 


-D(9(t))A x N = N(K(x) - N ) in fi, 
J^N(x,t) = 0 on 3D. 


( 29 ) 
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We can follow the same derivation as before, and discover that the value 0(t) of the fittest dispersal 
trait is now obtained through the time evolution constrained Hamilton-Jacobi equation 


'u t (8,t) - D(6)\V g u\ 2 = -H(6,N(;t )) 
max^ u(9 , t) = 0 = u(9(t),t). 


m 


(30) 


Here the effective fitness H(9,N(-,t)) is still given by the eigenvalue problem (16) with p = N. 

Note that, since derivatives vanish at a maximum point, we conclude that 

H(Q(t),N{;t))= 0 , 

and we also have N(x, t ) = A f[x, 9(t), iV( •, t)) , where N(x, t ) = lim e _^o Pe(x, t), 

We recall that, still formally, we can derive from (30) a canonical equation for the fittest trait 0{t) 
which takes the form, see mmmm, 


d * 


e(t) = -(- D 2 e u(m,t )) .VeH(9(t),AT(;t)). 


dt 


(31) 


We note that this equation also describes the fact that 8{t) will evolve towards smaller values of H 
and thus of the dispersal rate D. 

The main difficulties compared to the stationary case are to derive a priori estimates for p e (x,t ) 
analogous to those in Lemma [2.1| and to obtain gradient estimates on u £ . Since the Lipschitz regularity 
of u is optimal and only differentiability can be proved at the maximum point |S3], giving a meaning 
to (31) is also a challenge. 



Figure 1: Snapshots of the evolution of the dispersal for three different times according to the parabolic 
problem (28). Top figures depict the space repartition p{x) (fixed blue curve = K, increasing red = 
p). Bottom figures depict the trait distribution L 1 n(x, 9, t)dx. 
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Some numerics on the parabolic problem. Numerical simulations are presented in Figure [l] 
which illustrates the selection of lowest dispersal rate. Considering the parabolic problem (28), that 


means the convergence of the fittest trait 9(t) to the smaller values of the trait (with the coefficients 
below, this value is 0 = 0) as t —> oo. 

For this simulation we have chosen = (0,1) and the data 


D{9) = 1.5 6, K(x) = 1 + 20(1 - 4(s - .5) 2 )* 


and we have used, for the convenience of numerics, Dirichlet boundary conditions both in x and 9. 
The numerical scheme is the standard three points scheme in each direction, implicit in x and explicit 
in 9 because e = 10 -2 is small enough so as not to penalize the computational time. 

In Figure [lj we observe that the average trait, which initially is 9 ~ 0.7, decreases and gets close to 0 
at the third time displayed here, also the trait exhibits a concentrated pattern around its average. 
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